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1—5 ; ABSTRACT 

, If the high-energy emission from TeV blazars is produced by the Synchrotron Self- 

' Compton (SSC) mechanism, then simultaneous X-ray and gamma-ray observations of 

these objects are a powerful probe of the electron (and positron) populations respon- 
■ sible for this emission. Understanding the emitting particle distributions and their 

^ ' temporal evolution in turn allows us to probe physical conditions in the inner blazar 

0^ , jet and test, for example, various acceleration scenarios. Furthermore, by constraining 

(N . the SSC emission model parameters, such observations enable us to predict the intrin- 

sic (unabsorbed) gamma-ray energy spectra of these sources, a major uncertainty in 
current attempts to use gamma-ray observations to constrain the intensity of the Dif- 
fuse Extragalactic Background Radiation (DEBRA) at optical/infrared wavelengths. 
^-^ \ As a next step in testing the SSC model and as a demonstration of the potential power 

■~>^ I of coordinated X-ray and gamma-ray observations, we model in detail the X-ray and 

' gamma-ray light curves of the TcV blazar Mrk 501 during its April-May 1997 outburst 

^p^, with a time dependent SSC model. Extensive, quasi-simultaneous X-ray and gamma- 

Q ■ ray coverage exists for this period. We discuss and explore quantitatively several of 

\ the flare scenarios presented in the literature. We show that simple two-component 

C/3 , models (with a soft, steady X-ray component plus a variable SSC component) involv- 

ed ' ing substantial pre-acceleration of electrons to Lorentz factors on the order of 7niin = 

^ \ 10^ describe the data train surprisingly well. All considered models imply an emis- 

• 1^ . sion region that is strongly out of equipartition and low radiative efficiencies (ratio 

' between kinetic jet luminosity and comoving radiative luminosity) of 1 per-mill and 

^ \ less. Degeneracy in both, model variant and jet parameters, prevents us to use the 

- - I time resolved SSC calculations to substantially tighten the constrains on the amount 

of extragalactic gamma-ray extinction by the DEBRA in the relevant 0.5-50 microns 
wavelength range, compared to earlier work. 

Key words: galaxies: BL Lacertae objects: individual (Mrk 501) — galaxies: jets — 
X-rays: galaxies — gamma rays: theory 



(N 
O 



1 INTRODUCTION 

1.1 EGRET Blazar Observations 

The EGRET detector on board the Compton Gamma-Ray 
Observatory showed that many blazars are copious gamma- 
ray emitters (Hartman et al. 1999), their power at gamma- 
ray energies being comparable to (for low luminosity sources, 
i.e. BL Lac objects) or dominating by a wide margin (for 
high luminosity sources, i.e., FSRQs, Flat Spectrum Ra- 
dio Quasars, and OVVs, Optically Violantly Variables) the 
power emitted at longer wavelengths. The nonthermal ra- 



diation component probably originates from a population 
of relativistic particles embedded in the collimated outflow 
(jet) from a super-massive (lO'' up to several times lO^M©) 
black hole. The nonthermal continuum emission is com- 
monly explained with Synchrotron Compton (Ulrich et al. 
1997; Sikora & Madejski 2001) models: embedded in a jet 
which approaches the observer with relativistic velocity, a 
population of high energy electrons emits Synchrotron ra- 
diation at longer wavelengths and at shorter wavelengths. 
Inverse Compton (IC) radiation of high energy electrons off 
lower energy seed photons. The origin of the seed photons 
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is still uncertain (e.g. Blazejowski at al. 2000). The seed 
photon source could be "external" to the jet, e.g., radiation 
scattered and reprocessed by ambient matter in the Broad 
Line Region near the black hole, or infrared radiation emit- 
ted by dust in the inner nucleus of the host galaxy (External 
Compton models). Alternatively, the dominant seed pho- 
tons are synchrotron photons from the same electron pop- 
ulation responsible for the IC scattering (SSC, Synchrotron 
Self Compton models). In a generic source, both external 
and internal seed photons could be important in producing 
the observed spectrum. In the following we use the term 
Synchrotron Compton models if we do not want to specify 
the source of the seed photons. 

Alternative models, so-called "hadronic" models, invoke 
hadronic interactions of a highly rclativistic outflow which 
sweeps up ambient matter (Pohl & Schlickeiser 2000), inter- 
actions of high energy protons with gas clouds moving across 
the jet (Dar & Laor 1997), or, interactions of extremely high 
energy protons with ambient photons (Mannheim 1998), 
with the jet magnetic field (Aharonian 2000) , or with both 
(Miicke et al. 2002). If the reported fluxes of the diffuse in- 
frared background between 60 and 100 micron (Lagache et 
al. 1999, Finkbeiner et al. 2000) correctly describe the DE- 
BRA intensity in the far-infrared band, the "reconstructed" 
spectrum of Mrk 501, corrected for intergalactic absorption, 
may contain a sharp pile-up at and above 15 TeV. The latter 
cannot be explained by conventional Synchrotron Compton 
models. It has been argued that the presence of such a pile- 
up can be explained by bulk-motion comptonization (in the 
deep Klein-Nishina regime) of the ambient radiation by an 
ultra-relativistic conical cold outflow with a bulk Lorontz 
factor of ^10*^, while the remaining part of the spectrum 
could be explained by a conventional SSC model (Aharo- 
nian et al. 2002). 

All these models have some degree of success in ex- 
plaining the overall Spectral Energy Distribution (SED) of 
gamma-ray blazars. However, one can break much of the ap- 
parent degeneracy between these models by talcing advan- 
tage of the rapid, large-scale time variability these sources 
exhibit. Different models, for example, produce emission at 
a given frequency using particles of different energies, in- 
teraction cross-sections, and cooling times. The response of 
different models to changes in source conditions or the in- 
jection of fresh new particles is therefore different and in 
principle distinguishable - provided that one has sufficient 
time resolution to fully sample the flux variations and suffi- 
cient frequency coverage to constrain the different emission 
components that may be present. 

In view of this potential payoff, considerable effort has 
been dedicated to carrying out multi-wavelength observa- 
tions on powerful EGRET blazars like 3C 279 (Wehrle et al. 
1998). While the campaigns have lent considerable support 
to Synchrotron Compton models, the results of the cam- 
paigns wore not as conclusive as one might have hoped. The 
reasons for this are three-fold: 

(i) These blazars turned out to be highly variable on 
timcscalcs down to at least hours (Mattox ct al. 1997; Wag- 
ner et al. 1997). Even for the brightest objects, the instru- 
ment available for the gamma-ray observations, EGRET, 
simply did not have enough collection area to track all the 



gamma-ray flux variations, let alone provide high quality 
energy spectra. 

(ii) In typical models the electrons responsible for the 
GeV EGRET IC flux emit their synchrotron radiation at 
~UV energies. However, UV observations are difficult if not 
impossible because of atmospheric and galactic absorption. 
Thus the simultaneous observations that were made, e.g., 
at gammarray and X-ray energies, tracked radiation from 
electrons with very different energies and different cooling 
times and thus potentially different time histories and per- 
haps even emission regions. 

(iii) The observations showed that the gamma-ray emis- 
sion in several EGRET blazars is not consistent with the 
SSC model, the simplest version of Synchrotron Compton 
models (sec e.g. the comprehensive modeling of 3C 279 
broadband data described by Hartman et al. 2001). The ne- 
cessity to consider in External Compton models alternative 
seed photon fields substantially complicates the unambigu- 
ous interpretation of the data, especially since along our line 
of sight the beamed emission from the jet often dominates, 
making direct observation of these other photon fields diffi- 
cult. 



1.2 Potential of TeV Blazar Observations 

The second class of gamma-ray emitting blazars that 
EGRET discovered, the low power BL Lac objects like 
Mrk 421, were initially passed over as targets for extensive 
multi-wavelength campaigns since they were too weak in 
the EGRET band. The arrival of ground-based gamma-ray 
detectors like Whipple, HEGRA, and CAT with detection 
areas on the order of 10^ m^, however, now allows us to 
follow their gammarray fluxes on minute timescales (Gai- 
dos et al. 1996) and to routinely obtain detailed spectral 
information on timcscalcs down to one hour (Aharonian et 
al. 1999a). Besides their better accessibility at gamma-ray 
energies, these low power objects have several other impor- 
tant advantages. BL Lacs and their likely FR-I radio galaxy 
parent population appear to have underluminous accretion 
disks, i.e., "external" photon fields may not be important 
as seeds for IC scattering (Chiaberge et al. 1999). This to- 
gether with the fact that their time-averaged SEDs have 
successfully been described with one-component SSC mod- 
els, strongly suggest that SSC models which have much fewer 
free parameters than External Compton models indeed ap- 
ply. Also, perhaps because of the lower internal and external 
radiation fields and thus lower radiative losses (Ghisellini et 
al. 1998), the characteristic electron energies appear to be 
higher for the lower power objects, moving their synchrotron 
emission peak out of the UV, squarely into the X-ray range, 
where individual flares strongly dominate the overall lumi- 
nosity and can readily be observed with broad-band X-ray 
satellites like RXTE and BeppoSAX."" In the SSC model, 
the IC peak then moves from GeV to ~TeV energies. Thus, 
simultaneous X-ray and TeV gamma-ray observations follow 
the evolution of the electron population responsible for the 
bulk of the source luminosity, and the observations are well- 
matched in the sense that they track the emission from the 
same electrons, providing tight constraints on the electron 
distribution and its time evolution. SSC models that apply 
to these objects are therefore testable, especially with the 
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next generation of X-ray and gannna-ray detectors coming 
on line in the next few years. 

Proving whether an SSC model works or not has a po- 
tentially largo payoff. If the model does not work, then we 
must significantly revise our understanding of the physical 
conditions and processes in these objects. If it does work, 
then we can use it for example to probe the acceleration 
processes at work in the innermost region of the jet. We can 
also use it to constrain the amount of extragalactic gamma- 
ray extinction due to pair production processes on the difi'use 
optical/infrared background 7TeV + 7ir,o — » e"*" e~ (Gould 
& Schreder 1966; Stecker, De Jager & Salamon 1992), by 
comparing the predicted intrinsic TcV gamma-ray energy 
spectrum with the observed one (Coppi & Aliaronian 1999, 
Krawczynski et al. 2000, called "Paper I" in the following"). 

Based on these considerations, the brightest TeV 
blazars, Mrk 421 {z = 0.031) and Mrk 501 {z = 0.034), have 
been the subject of increasingly intensive observing cam- 
paigns. This has led to the discovery of pronounced TeV 
gamma-ray / X-ray flux correlations for Mrk 421 (Buckley 
et al. 1996, Takahashi et al. 1996, Maraschi et al. 1999, Taka- 
hashi et al. 2000) and Mrk 501 (Pian et al. 1998, Djannati- 
Atai et al. 1999, Paper I, Sambruna et al. 2000, see also 
Fig. 1 of this paper). 

1.3 Goal of this Paper and Relation to Previous 
Work 

The goal of this paper is to extend the analysis of Paper I, 
which was a first joint analysis of an unprecedented set of 
X-ray/TeV monitoring data taken during the 1997 flare of 
Mrk 501. Using RXTE (X-ray) and HEGRA (TeV) observa- 
tions that were simultaneous to within a few hours (i.e. less 
than the ~12 hour characteristic variability time scale of the 
source) , we showed that the X-ray flux of the source, partic- 
ularly above 10 kcV, was strongly correlated with the TeV 
flux, in accord with Synchrotron Compton models. More- 
over, we found that for that two month period the data were 
consistent with the quadratic relation expected between the 
X-ray and TeV fluxes in a simple SSC model, although a 
linear relation between X-ray and TeV flux plus a constant 
base X-ray flux level also described the data satisfactorily. 
We then used an one-zone, steady state SSC model to fit the 
X-ray /TeV energy spectra for several days in order to see if 
the model could explain the data and to make a first guess at 
the SSC source parameters. We found that the data could 
be fairly well-described by a reasonable sequence of SSC 
models. By assuming that we were indeed seeing SSC emis- 
sion and by taking (at the time) extreme values of the SSC 
model parameters (e.g., jet Doppler factors ~100), we then 
placed constraints on the maximum amount of extragalac- 
tic gamma-ray absorption present in the observed spectrum 
(Paper I, Fig. 10). We will give an updated very detailed 
discussion of the implications of SSC and External Comp- 
ton models on the intensity of the DEBRA in a companion 
paper (Coppi et al. 2001). 

In this paper, we attempt to quantify how well SSC 
model predictions match observations by taking the 1997 
data set and fitting the full observed spectral sequence us- 
ing a time dependent SSC code (Coppi 1992). The code 
accurately models the temporal evolution of the energy dis- 
tribution function of a population of relativistic electrons 



due to acceleration processes and radiative and adiabatic 
energy losses. However, it is a one-zone code that assumes 
homogeneous and isotropic particle and pitch angles dis- 
tributions in the jet rest frame. Leaving aside plasma and 
magnotohydrodynamics issues, such assumptions are clearly 
an oversimplification given the inhomogeneous structure of 
jets, especially as observed on VLBI radio scales where sev- 
eral distinct "blobs" (emission regions) may be active at any 
given time. However, during 1997 the Mrk 501 emission was 
strongly dominated by individual flares during which the X- 
ray flux increased by up to 5 times and the TeV gamma-ray 
fluxes by up to 30 times. It is highly probable that these indi- 
vidual flares were produced by single jet regions, rather than 
being the superposition of several, causally not connected 
events. In addition, the physical conditions in all the emis- 
sion regions seemed to be very similar: the X-ray and TeV 
gamma-ray fluxes were well correlated during more than two 
months (sec Paper I, and this paper Fig. 1), and the TeV 
energy spectrum stayed remarkably stable during more than 
6 months (Aharonian 1999arc). As a final justification of our 
approach, the analysis presented in Paper I showed that the 
X-ray and TeV gamma-ray fluxes typically varied on time 
scales of ~1 day with shortest flux rise and decay times on 
the order of half a day. Kataoka et al. (2001) and Tariihata 
et al. (2001) analyzed RXTE and ASCA data taken dur- 
ing the years 1997-2000 and find, in accord with our results, 
that Mrk 501 has a low duty cycle for flares on time scales 
of a few hours and shorter. Thus, the sampling of the data 
with 2 X-ray observations and several TeV gamma-ray ob- 
servations per day was probably sufficient for giving a rough 
picture of how X-ray and TeV gammar-ray fluxes evolved in 
time. 



Observational signatures of Synchrotron Compton mod- 
els have been described by various authors (see e.g. the ref- 
erences given in Table 1). In the following we show for the 
first time an attempt to fit a prolonged sequence of X-ray 
and TeV gamma-ray data with a time dependent SSC code. 
This approach makes it possible to use the full information 
encoded in the correlated flux variability at different wave- 
lengths. In contrast to parametric SSC fits (see e.g. Paper 
I, Tavecchio et al. 2001) the method uses a self-consistently 
evolved electron population which assures that the assumed 
electron energy spectrum is physically realizable from an ini- 
tial acceleration spectrum (see also the discussion by Mas- 
tichiadis & Kirk 1997). Wo think that the approach of us- 
ing a time resolved analysis to break model degeneracies 
will become increasingly more powerful and important as 
the sensitivity and energy coverage of X-ray and gamma- 
ray instruments continue to improve. Note that a thorough 
understanding of the SSC model is also a necessary prereq- 
uisite for the evaluation of External Compton models which 
always include a SSC component. 



The rest of this paper is structured as follows. In Sect. 
2 we introduce the data set and show an updated version of 
the X-Ray/TeV gamma-ray flux correlation. In Sect. 3 we 
describe the model calculations and in Sect. 4 the time de- 
pendent model fits. Finally, we discuss the results in Sect. 5. 
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2 THE DATA SET 

During 1997 the BL Lac object Mrk 501 went into a re- 
markable state of continuous strong flaring activity and the 
source was intensively monitored at X-rays and TeV gamma- 
rays. During April and May, 1997 the source was regularly 
observed with the RXTE X-ray satellite, with typically two 
pointed observations per day of ~20 min duration (Paper I). 
Each pointing resulted in a high accuracy measurement of 
the 3-25 keV X-ray flux and photon index with an accuracy 
which was only limited by systematic effects. The curva- 
ture of the X-ray spectrum could be assessed for a couple of 
pointings with relatively high X-ray fluxes and long integra- 
tion times. On three days (April 7, 13th, and 16) the source 
was also scrutinized with the BeppoSAX X-ray telescopes, 
revealing the X-ray energy spectrum of the source over the 
broad energy range from 0.1 keV to ~ 200 keV (Plan et al. 
1998). 

In Paper I we studied the correlation of the X-ray 
fluxes with the TeV gamma-ray fluxes as measured with 
the HEGRA Cherenkov telescope system (Aharonian et al. 
1999a). For the present study we complemented the data 
set with the TcV fluxes from the HEGRA CTl (Aharo- 
nian et al. 1999c), Whipple (Quinn ct al. 1999), and CAT 
(Djannati-Atai et al. 1999) telescopes. In Paper I we found 
a very tight correlation between the 25 keV and 2 TeV 
fluxes. The flux variability amplitude was approximately 3 
times larger at TeV than at X-ray energies, being consis- 
tent with a quadratic relationship. An updated version of 
the X-ray/TeV gamma-ray flux correlation is shown in Fig. 
1. The additional X-ray /TeV gammar-ray flux pairs confirm 
the previous finding of a clear fiux correlation. However, the 
quality of the correlation still does not allow us to differenti- 
ate between a quadratic X-ray/TeV gamma-ray relationship 
and a linear one with a non-zero X-ray fiux offset. 

The X-ray as well as TeV gamma-ray data are plagued 
by systematic errors. In the case of the BeppoSAX data the 
spectral index below 1 keV is not well determined due to un- 
certainties in the neutral hydrogen column density. Above 
50 keV the scatter of the data points increases more than the 
statistical errors, indicating systematic uncertainties in the 
detector response and/or the background subtraction proce- 
dure. Some TeV gamma-ray points taken at approximately 
the same time with different experiments deviate by more 
than 3 a statistical error from each other, indicative either 
for very fast source variability, or, for errors due to unstable 
atmospheric conditions. While the majority of nearly coinci- 
dent measurements shows good agreement between different 
TeV telescopes, the occurrence of some exceptions makes it 
difficult to decide between models if their predictions dif- 
fer only for one or two days. This caveat will be discussed 
further below, when we compare the models with the data. 

The spectral variability at TeV energies has been a mat- 
ter of debate: the HEGRA group did not detect spectral 
changes with an accuracy (1-5 TeV photon index) of ~ 0.2 
and 0.05 for diurnal and flux selected mean energy spectra, 
respectively. The CAT group reported the statistically sig- 
nificant detection of a hardness intensity correlation based 
on the F(> 900GeV)/F(> 450 GeV) hardness ratio, corre- 
sponding to a ~0.25 change in photon index. The two data 
sets overlapped only partially in time: the HEGRA group 
did i.e. not take data on April 16, 1997, which is the most 
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Figure 1. Correlation between X-ray (RXTE) and TeV gamma- 
ray fluxes. Tlie gamma-ray fluxes are from CAT (squares), 
HEGRA CT System (solid points), HEGRA CT 1 (asterisks), 
and Whipple (open circles) . Only observation pairs with less than 
6 hrs time delay have been used. 



important day in the CAT analysis. Koriopelko ct al. (1999) 
noted that the stability of the TeV energy spectra, evident 
in the HEGRA data, might be used to infer constraints on 
the intensity of the DEBRA. In the plots shown below, we 
cross-calibrated the BeppoSAX data relative to RXTE mea- 
surements taken at approximately the same time. 

Compared to the results shown by Plan et al. (1998) we 
reduce the normalization of the BeppoSAX PDS data by up 
to 35% which eliminates the discontinuity of the joint Bep- 
poSAX MECS, LEGS, and PDS energy spectra at ~15 keV 
(between the energy coverage of the LEGS and PDS instru- 
ments) and is then consistent with the spectral shapes si- 
multaneously measured from 3 keV to 25 keV with RXTE. 
Wc also cross-calibrated the CAT, HEGRA CTl, and Whip- 
ple gamma-ray fluxes relative to the ones measured by the 
HEGRA CT System. Although we obtained a list of CAT 
fluxes as function of the integer MJDs of the observations, 
the fractional MJDs of the CAT observations are not known 
to us. In the following we centered the CAT observations at 
12 am UTC. 



3 MODELING 

3.1 Time Dependent SSC Code 

The SSC code (Coppi & Blandford 1990; Coppi 1992) as- 
sumes a spherical emission region of radius R which is flUed 
with an isotropic electron population and a randomly ori- 
ented magnetic fleld B and which approaches the observer 
relativistically. The motion of the jet toward the observer 
can be characterized with the jet Doppler factor, defined by 



= r(i-/?cos(e)). 



(1) 
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with r the bulk Lorentz factor of the emission plasma, and 
P its bulk velocity in units of the speed of light, and 9 is 
the angle between jet axis and the line of sight in the ob- 
server frame. The TeV gamma-ray flux variability on time 
scale ATobs ~ 12 hr (Aharonian ct al. 1999a) together with 
causality arguments set an upper limit on the radius of the 
emission volume: 



(2) 



If the jet moves along a curved path more rapid flares could 
result from a change of the Doppler factor as the jet's radi- 
ation beam sweeps across the observer. 

The kinetic equations, discrctizcd in energy, take fully 
into account the non-continuous character of IC processes in 
the Klein-Nishina regime, and are evolved in time with a two 
step implicit scheme treating first the photon distribution 
and subsequently the electron distribution. The length of 
time steps is chosen such that the number of photons and 
particles per energy bin changes per step by less than 20%. 
The kinetic equation for the photon density (per unit volume 
and energy) n-y reads: 



'dt 



= q-i - Pi - 



R{1 + k) 



(3) 



where q-y de and p-, de are the rate of photons being produced 
into and out of the cncrgj' interval [e, e + de] due to electron- 
magnetic field, electron-photon and 2-photon interactions. 
The last term of the right hand side represents photons 
which escape from the emission region. The factor ci?~^in 
the last term assures that the photon density approaches 
steady state values only with a rise/decay constant longer 
than the light crossing time. The factor (1-|-k(7)) parameter- 
izes the modification of the photon escape time by Compton 
processes (Coppi 1992); however, for all the models discussed 
in the following, we have always k -C 1. The Klein-Nishina 
effect decisively influences the resulting gamma-ray energy 
spectrum and proper modeling is imperative. 

The kinetic equation of the electron (and possibly 
positron) density rie reads: 



One ^ dr. , 
~dt ^ 97 



Tie 



(4) 



with Qe(7, t) from Eq. (5), 7cont gives the decrease of an elec- 
tron's Lorentz factor per unit time due to continuous energy 
losses, and qe d'y and pe d'y are the rate of particles being pro- 
duced or scattered into and out of the Lorentz factor interval 
[7, 7 + dj] due to non-continuous energy loss processes, re- 
spectively. The last term of the right hand side represents 
an energy independent escape probability of electrons from 
the emission region. 

To first order approximation our code takes the non- 
vanishing source extension into account through the last 
term in Eq. 3. As a consequence, the code is able to describe 
flux variations oven on time scales on the order of _R/c in a 
qualitatively correct way. We limit ourselves in this paper to 
describe the time variable emission component with a one- 
zone SSC model. A one-zone model is able to approximate 
multi-zone models as long as the spatial gradients of the 
magnetic field and the non-linear components in the prop- 
erly modified kinetic equations (3) and (4) are small. Our 
code can i.e. mimic "linearized inhomogeneous models" as 
discussed by Kirk, Rieger & Mastichiadis (1998) and Chi- 



aberge & Ghisellini (1999). While External Compton mod- 
els can be dominantly linear, the SSC model is inherently 
non-linear: the synchrotron component directly follows the 
evolution of the electron population, but the IC component 
results from the interaction of the electron population with 
the self-produced synchrotron photons. Since electrons and 
synchrotron photons traverse the emission region on a time 
scale of R/c, one expects that the IC component lags the 
synchrotron component by approximately one light crossing 
time (Coppi & Aharonian 1999). This is the most drastic 
time lag effect expected in the SSC model. For Mrk 501 how- 
ever no such time lag has been observed so far, the upper 
limit being about 12 hrs (Aharonian ot al. 1999a; Aharo- 
nian et al. 1999c; Paper I; Sambruna et al. 2000). As long as 
instrumental resolutions do not permit to resolve this time 
lag, we think it is safe to use only one component to describe 
the time variable emission. 

We flt the full two months data train using a single 
emission volume. As we will point out in the discussion, it 
might be that individual flares (of durations on the order 
of '^1 day) are produced by independent emission regions. 
Upon flaring, a region would expand adiabatically, and thus 
fade away quickly. Even in this case, our model should give 
reasonable results for two reasons: (i) as it turns out the 
best fitting models have particle escape times on the order 
of the fiux variability time scale; (ii) the tight X-ray/TeV 
gamma-ray flux correlation argues for a very similar size of 
the emission regions. As a consequence, each flare is pro- 
duced by freshly accelerated electron populations and mod- 
eling the flares with one emission region gives similar results 
as using several disjunct emission regions. 

3.2 Treatment of Particle Acceleration 

Given the sparse observational sampling of our data set in 
time and wavelength, we did not embark on modeling the 
acceleration process in detail but used instead an "exter- 
nal" acceleration function. We parameterize the production 
rate of freshly accelerated particles as function of electron 
Lorentz factor 7, spectral index of particle acceleration p, 
normalization Qo{t), minimum Lorentz factor 7min, and high 
energy cut-off 7max(i) as follows: 

Qe(7, t) = Qo{t) exp (-7/7max(t)) 6(7 - 7mi„) (5) 

with e(a;) = for x < and e(x) = 1 for a; > 0. We 
use the canonical value of p = 2 expected for diffusive par- 
ticle acceleration at strong shocks (Bell 1978; Blandford & 
Ostriker 1978) and do not consider the ramifications arising 
from the non-linear modification of the shock structure due 
to the backreaction of accelerated particles (Bell 1987) and 
mildly or ultra-relativistic shock velocities (see the recent 
review by Kirk & Duffy 1999). 

The low-energy cutoff in the spectrum of accelerated 
electrons 7min is a critical model parameter. If the radiative 
cooling time of electrons with Lorentz factor 7min is shorter 
than all the other characteristic time scales of the system, 
the main break of the electron spectrum occurs at 7min. 
Thus, at high enough values (~ 10^), 7min determines the en- 
ergies at which the synchrotron and IC SEDs peak. On the- 
oretical grounds one expects much lower values of between 
1 and the proton to electron mass ratio mp/rrie = 1836 
(Hoshino et al. 1992; Levinson 1996; McClements 1997). We 
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will use in the following a relatively low value of 7min = 
1000 as the fiducial value and will discuss higher values at 
several points. 

We characterize the acceleration luminosity le by the 
pair-compactness parameter (Coppi 1992): 

3.3 Treatment of Extragalactic Extinction 

The TeV gamma-ray spectra are expected to be modified 
by extragalactic extinction due to pair production processes 
of the TeV gamma-rays with photons of the DEBRA. The 
uncertain DEBRA intensity in the relevant 0.5-50 microns 
wavelength range introduces a major uncertainty in the 
modeling of the source. While earlier estimates of the DE- 
BRA level predicted negligible extinction at gamma-ray en- 
ergies below ~1 TeV, more recent observational and theo- 
retical efforts suggest that this might not be true (Primack 
et al. 2001). We think that model estimates of the DEBRA 
still have not reached the reliability that we should limit 
our computations to a specific DEBRA model. Rather we 
will treat the modification of the TeV fiux level and energy 
spectrum as not fully constrained. Clearly, the DEBRA ex- 
tinction does not modify the relative TeV gamma^ray flux 
variations and we use the information encoded in the relative 
fiux changes by fitting the TeV gamma-ray fiuxos subject to 
a common constant scaling factor ^. At 2 TeV one expects 
a ^-value of between 0.2 and 1. We have varied the ^-values 
in this range, and the qualitative conclusions presented bo- 
low are robust and do not depend on the exact value of ^. 
For the detailed fits presented below, wo will take ^ ~ 0.5, 
a value which seems to be preferred by recent observations 
and theoretical modeling. 

Due to rather large statistical errors on diurnal gamma- 
ray photon indices we did not attempt to fit the variations 
of the TeV energy spectra with the SSC code. We did check 
that the modeled TeV gammar-ray energy spectra are con- 
sistent with the observed one, taking into account that ex- 
tragalactic extinction only steepens the gamma-ray energy 
spectra. 

3.4 Fitting Procedure 

The free parameters of our model are the radius of the emis- 
sion volume R, the jet Doppler factor 5j , the mean magnetic 
field B, the escape time of relativistic electrons from the 
emission region fosc , the normalization of the electron accel- 
eration rate Qo, and the minimum and ma^dmum Lorentz 
factors of accelerated particles 7min and 7max. We fit the 
April and May, 1997 RXTE 10 keV fiuxcs and 3-25 koV 
photon indices and the 2 TeV fiuxes derived from CAT, 
HEGRA, and Whipple measurements. Given a hypothesis 
of what causes the flaring activity (a variable Qo{t), 7max(t), 
and/or 5j{t)) we fit the data in 2 steps: 

(i) For a set of parameters {R, Sj, B, 7min, tcsc) we de- 
termine the simplest possible function Qo{t) and 7max(t), or 
for some models Sj{t), which fit the X-ray flux amplitudes. 
Hereby, "the simplest possible functions" means that given 
the X-ray flux measurements at times ti (in the jet frame) 



we use simple prescriptions to determine Qo, 7max, or 5j at 
the times t[ — ti — At, and compute intermediate values 
by simple interpolation. The parameter At is the time by 
which the photon density in the emission region reacts to 
changes of the electron spectrum. In general, the optimal 
delay At depends on the time scale on which electrons cool 
and escape from the emission region. Due to the observa^ 
tional constraint on the time lag between low (3 keV) and 
high energy X-rays (30 keV) to be shorter than '^10 hrs 
(Paper I), the delay of all our models is dominated by the 
light crossing time and the treatment of Eq. (3) results in 
an optimal value of At ~ 2.5 iic^^. For models in which 
ordy 5j(t) produces the time variability we use At = for 
obvious reasons. 

We determine the values Qo{ti), 7max(ti), and 5j(ti) iter- 
atively by making a first guess, computing the SSC model, 
and adjusting the values until the X-ray fluxes are described 
satisfactorily. Usually, between 2 and 5 iterations are needed. 
The reader should keep in mind that the true time history 
of the acceleration process could be more complex. 

(ii) We vary the parameters {R, Jj, B, 5^, tesc) to obtain 
the best fit to the observed X-ray photon indices and TeV 
gamma-ray flux amplitudes. The quality of the fits is charac- 
terized by x^-values, computed for the X-ray photon indices 
and for the TeV gamma-ray fiuxes. We exclude the first 2 
days of each observation period from entering the x^-values, 
since the results strongly depend on the unknown behav- 
ior of the source before the observations commenced. Note 
that the reduced x^-values of the X-ray fiuxes and photon 
indices and of the TeV gamma^ray fluxes exceed 1 by a wide 
margin, showing that the experimental statistical errors axe 
smaller than the accuracy of our models and/or that the ex- 
pcrimcutal systematic errors (often only poorly determined) 
are non-negligible. 

Computing the -values, we scale all TeV gamma-ray 
fiuxes by a common factor ^ with 0.2 < ^ < 1. For each 
model we state the ^-value which we used as well as q. The 
latter value is the difference of the observed and the mod- 
eled mean spectral index. We interpret this difference as due 
to extragalactic extinction. Due to the strong dependence 
of the IC luminosity on the radius of the omission volume, 
the modeling does not give any constraints on ^. As will 
be discussed in Sect. 5, the uncertainty in the SSC model 
parameters do not allow us to constrain rj either. 

After obtaining in this way the best fits we test the predicted 
SEDs (Spectral Energy Distributions) for consistency with 

the broadband X-ray spectra from the BeppoSAX observa- 
tions, and an MeV upper limit from EGRET (Catanese et 
al. 1997). 



4 RESULTS OF THE TIME DEPENDENT 
MODELING 

SSC blazar models have extensively been discussed in the 
literature. The models can roughly be classified according 
to 2 criteria (see Table 1): 

(i) According to what produces the observed gamma-ray 
flares: basically, almost every parameter of the SSC model 
has been invoked by at least one group to account for the 
blazar flaring activity. 
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(ii) According to the mechanism that determines the en- 
ergies at which the synchrotron and IC SEDs peak. The SED 
peak energies are either determined by the minimum Lorentz 
factor 7min of accelerated particles, or, by the balance be- 
tween radiative cooling times and the shorter of particle es- 
cape time and the characteristic duration of individual flares 
(sometimes referred to as injection time scale, or, dynamical 
time scale of the jet). 

The modeling of the full data train is computationally very 
intensive and we therefore focused on exploring only the 

models which seemed most promising to us. While the time 
resolved analysis clearly rules out some models, it gives fits of 
very similar quality for others. Our difflculties to distinguish 
between models mainly derive from two facts: 

(i) From the limitations of the data set, namely sparse 
observational sampling in time and energy, and systematic 

errors on X-ray energy spectra and TcV gamma-ray fluxes. 

(ii) From the unknown modiflcation of the TeV gamma- 
ray energy spectra by extragalactic extinction. 

Keeping these limitations in mind, we discuss the fit results 
with a focus on pointing out which models are capable of 
correctly describing the qualitative behavior of the X-ray 

and TcV gamma-ray radiation. Table 2 lists the model pa- 
rameters of the SSC models shown in the figures. 

4.1 One-Component Models 

4.. 1.1 Time variability through Qo{t) 

Wo first consider time variability through a varying rate of 
accelerated particles. If Qo(i) varies, the SSC mechanism 
automatically produces a more than linear increase of TeV 
flux as function of X-ray flux. Fig. 2 shows the observed and 
modeled X-ray and gamma-ray flux amplitudes and pho- 
ton indices. In Paper I, wo derived a lower limit on the 
Doppler factor of 6.3. In most of the following models we 
use a rather high Doppler factor of Sj = 45, for two reasons: 
(i) a high Doppler factor allowed us to fit the data with a 
wide range of magnetic field values; for lower Doppler fac- 
tors, weak magnetic fields result in a strong overproduction 
of TeV gamma-rays (see the related discussion by Krawczyn- 
ski et al. 2001); (ii) non- negligible extragalactic extinction of 
TeV gamma-rays seems highly probable; the high Doppler 
factor results in predicted TeV energy spectra that agree 
with the observed ones after correcting for extragalactic ex- 
tinction (which steepens the TeV gamma-ray energy spec- 
tra) . Ideally, we would like to use the modeling to determine 
the amount of extinction. As we will discuss further below, 
for the time being we can not do this due to parameter de- 
generacies. 

Whenever the size of the emission region satisfies 
Eq. (2), and either rapid particle losses (tesc not much larger 
than Rc~^) or a sufficiently large magnetic field allow flares 
to decay rapidly enough, the X-ray amplitude can be de- 
scribed to arbitrary precision. This also applies for all the 
models described in the following. Thus, we subsequently 
focus on the X-ray photon indices and the TeV gamma-ray 
fluxes for measuring the quality of a fit. 

Although the model describes the TeV flux levels, it fails 
to reproduce the range of observed X-ray photon indices. 
The stability of the X-ray spectrum is a solid property of this 



model. The spectral variability shown in Fig. 2 is already the 
result of a fine tuning between the parameters B, tcsc, and 
7max- Most realizations of this model result in substantially 
less spectral variability. For the model parameters of Fig. 2 
the location of the break in the synchrotron SED is given 
by the competition of the escape of electrons on time scale 
tesc and their radiative cooling through synchrotron and IC 
emission on time scale trad- High energy electrons with trad ~ 
Rc~^ cool almost instantly. Low energy electrons with trad ^ 
tcsc do not have time to cool before they escape the emission 
region. The result is that the spectrum does only vary over 
a rather small region where trad ~ iesc- The "smearing" of 
the break in the electron spectrum duo to the width of the 
synchrotron emissivity results in a rather stable break of 
the synchrotron spectrum. Qualitatively, we see a similar 
behavior also for larger values of tesc ■ The reason is that the 
rise and decay times of the flares are of the same order of 
magnitude as the time between flares. Thus, all observed 
electron energy spectra radiatively cool to approximately 
the same degree with a break at approximately the same 
Lorentz factor. 

Inspection of the fitted and observed TeV fluxes shows 
a strong discrepancy at MJD 50544. A TcV data point that 
suggests a very low flux is bracketed by X-ray observations 
of relatively high fluxes. This model as well as the mod- 
els described below fail to describe this exceptional anti- 
correlation. The low TeV flux may either be indicative of a 
short period of low X-ray and TeV gamma-ray activity be- 
tween the two X-ray observations, or of an underestimated 
TeV flux due to instrumental or atmospheric irregularities. 
Note how little the TeV photon index changes in this model. 

4.1.2 Time variability through 7max(t) 

As a second model we tested flares caused exclusively by a 
variation of the high energy cutoff 7max of accelerated parti- 
cles. While 7max may depend on the details of the magnetic 
field structure in the surrounding of a particle accelerating 
shock, other parameters influencing mainly the acceleration 
of lower energy electrons could remain constant. Historically, 
such models were motivated by observation of the blazar 
Mrk 421 which showed dramatic X-ray and TeV gamma- 
ray flux variability accompanied by only minor optical flux 
variability. 

Varying 7max alone we did not achieve a satisfactory fit 
to the data. A typical result is shown in Fig. 3. As before the 
X-ray fluxes can be described to arbitrary precision. While 
the model has no difficulty in producing the observed range 
of X-ray photon indices, it fails to describe the X-Ray / TeV 
gamma-ray flux correlation: the predicted TeV gamma-ray 
fluxes hardly vary at all. Combinations of Doppler factor and 
magnetic field where the TeV fiux changes more strongly 
than the X-ray flux result in steeper than observed TeV 
energy spectra. The time variation of 7max causes largo flux 
and spectral variability only at energies ^10 TeV where 
the inherent TcV energy spectrum is extremely soft (photon 
index ^ 2.5). Extragalactic extinction can not remedy this 
shortcoming since it is believed to steepen and not to soften 
the TeV energy spectra. 

Tavecchio et al. (2001) studied parametric SSC model 
fits to Mrk 501 snapshot data and concluded that the max- 
imum Lorentz factor of accelerated particles is mainly re- 
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Figure 2. X-ray and TeV gamma-ray data (data points) from 
April, 1997 (left column) and May, 1997 (right column) with SSC 
model fit (lines). Here, model flares are caused by a time depen- 
dent Qo{t). The panels show from top to bottom: (i) the logarithm 

of the injection compactness le (sec Eq. 6), (ii) the logarithm 
of the 10 kcV X-ray energy fiux (CGS units), (iii) the X-ray 3- 
25 keV photon index, (iv) the logarithm of the 2 TeV energy flux 
(CGS units), and (v) the 1-5 TeV photon index. The gamma- 
ray fluxes are from CAT (squares), HEGRA CT System (solid 
points), HEGRA CT 1 (asterisks), and Whipple (open circles). 
TeV photon indices have only been published by the HEGRA 
group. The value for April 16 has been inferred from the energy 
spectrum published by Djannati-Atai et al. (1999). The vertical 
dotted lines show the days with BcppoSAX observations (April 7, 
11, and 16 which will be discussed in more detail further below). 
The model parameters arc: 5j = 45, R = 1.1 X 10^'' cm, B = 



0.014 G, tesc 

0.5, T) = 0.2. 



10 i?c-i, 



In 



10^ 7max = 2.5 X 10', e 



sponsible for the flaring activity. However, detailed inspec- 
tion of their fit parameters shows that they described the 
data by varying the break of the electron spectrum rather 
than the high energy cutoff. Furthermore, the fits involve 
a substantial variation of the normalization of the electron 
spectrum (i.e. the acceleration rate). 



4.1.3 Time variability through Qo{t) and 7max(t) 

Models in which both, Qo and 7max, change with time 
have been invoked to account for the secular changes of the 
Mrk 501 X-ray (Plan et al. 1998; Sambruna ct al. 2000) and 
TeV gamma-ray (Aharonian et al. 2001a) energy spectra. In 
models of diffusive electron acceleration at strong shocks the 
electron acceleration rate Qo is determined by the rate with 
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Figure 3. Same as in Fig. 2, but with flaring activity through a 
time dependent maximum Lorentz factor of accelerated particles 
7inax(t). The upper panel shows here 7max(t). The model param- 
eters are: 5j = 45, i? = 1.5x 10^® ci 
KC-I, 7niin = 103, ^ = 0.5, jj = 0.2. 



which particles are "injected" into the acceleration process. 
The high energy cutoff of accelerated electrons 7max is deter- 
mined by the competition between electron energy gains and 
energy losses. Changing plasma properties most probably af- 
fects both, the injection rate and the high energy cutoff. Due 
to the uncertain nature of the particle injection mechanism 
we choose a simple parametric description to describe the 
correlation between Qo and 7max: 



7ii 



Qo" 



(7) 



and treat the exponent a as an additional free parameter 
of the fit. The x^-value of the X-ray photon indices show a 
pronounced minimum for a value of a = 2, and Fig. 4 shows 
a SSC fit to the data. The model describes the RXTE and 
gamma-ray data rather satisfactorily. 

During the first 3 days of the April campaign the model 
X-ray indices are by ~0.1 harder than the observed ones, a 
discrepancy which is shared also by all subsequent models. 
Compared to other days of similar X-ray flux levels, the X- 
ray spectrum of the first three days with RXTE coverage was 
very soft, indicating that the source properties did evolve 
during the 2 months campaign. 

Although only 3 BeppoSAX observations were per- 
formed during 1997, the data is very constraining since it 
covers the broad energy range from 0.1 keV to ~200 keV. 
The long BeppoSAX pointings of ~12 hrs duration brack- 
eted the RXTE and TeV gammarray observations. The com- 
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Figure 4. Same as in Fig. 3, but with flaring activity through 
time dependent Qo (*) and 7max (t) with 7max oc Qo ^ • The upper 
panel shows here Qo(*)- The model parameters are: <5j = 45, R = 
3.2 X IQis cm, B = 0.035 G, tesc = 3 Re"!, 7min = 10^, € = 
0.5, ?? = 0.1. 

parison of the modeled and observed broadband energy spec- 
tra is shown in Fig. 5. For all three days with BeppoSAX 
data the model fails to predict to observed X-ray energy 
spectra at and above 50 keV: the modeled energy spectra are 
all too soft. The high energy spectrum however is a rather 
solid prediction of this model: the cutoff at relatively low 
energies is needed to correctly account for the spectral vari- 
ability observed with RXTE. 

In Fig. 5 also the modeled and observed TeV energy 
spectra are shown. We used for all three days the time av- 
eraged 1997 HEGRA TeV gammarray energy spectrum nor- 
malized at 2 TeV to the mean flux measured on that day 
with all operational TeV telescopes. The use of the time 
averaged energy spectrum is justified by the fact that all 
the HEGRA data taken during 1997 are consistent with the 
shape of the time averaged energy spectrum (see however 
Djannati-Atai et al. (1999)). For illustrative purposes the 
dashed linos in Fig. 5 show the predicted TeV gamma-ray 
energy spectra modified by extragalactic absorption. We will 
further discuss the agreement between observed and mod- 
eled SEDs below for the models which give a better overall 
fit to the X-ray data. 

4.2 Two-Component Models 

The region near the presumed central black hole as well 
as various regions along the jet might emit X-rays in 




log(E/eV) 



Figure 5. For the model of Fig. 4 the upper panel shows the 
electron energy spectra {E^ dN/ dE, energy E in the jet frame) 
responsible for the observations of April 7 (solid line), April 11 

(dashed line) and April 16 (dotted line). The lower three pan- 
els compare the observed (points) with the modeled (solid line) 
SEDs (energy in observer's frame). For illustrative purposes the 
dashed line shows the TeV gamma-ray energy spectra modified 
by extragalactic absorption as predicted by the DEBRA model 
"LCDM, Salpeter Stellar Initial Mass Function" of Primack et al. 
(2001). For April 7, the model under-predicts the TeV flux, and 
the dashed-dotted line shows the same absorbed gamma-ray spec- 
trum, but normalized to the flux at 2 TeV to facilitate comparison 
of the spectral shapes. Asterisks show the BeppoSAX data nor- 
malized to the flux measured with RXTE. The solid points show 
the shape of the HEGRA 1997 time averaged Mrk 501 energy 
spectrum (Aharonian ct al. 1999b; Aharonian et al. 2001b) nor- 
malized to the mean TeV gamma-ray flux measured on each day 
by GAT, HEGRA, and Whipple. The 2a upper limit at 100 MeV 
has been derived from EGRET observations between April 9th 
and April 15th, 1997 under the assumption of a constant emis- 
sion level (Catanese et al. 1997). 



Bremsstrahlung, IC, or Synchrotron processes without pro- 
ducing a comparable luminosity at gamma-rays (see also 
Blazejowski et al. (2000), Bicknell, Wagner & Groves 
(2001)). Thus it is conceivable, even probable, that the X- 
ray emission from Mrk 501 is "contaminated" by an emission 
component which varies on longer time scales than the TeV 
gamma-ray radiation. While all the one-component mod- 
els described in the previous section failed to fully describe 
the data, we find that the addition of a quasi-stationary X- 
ray component substantially improves the situation. Vary- 
ing contributions of the quasi-stationary soft and the time 
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Figure 6. Fit of a SSC model with two emission components: (i) 
a quasi-stationary X-ray component, and, (ii) a time variable X- 

ray/TcV gamma-ray component, flares produced through Qo{t)- 
Data and units are the same as in Fig. 2. The model parameters 



are: 5; = 45, i? = 3.4 X 10^'' cm, B = 0.014 G, U. 
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Figure 7. Same as in Fig. 5 but for the two-component model 
with time variable Qoj(t) shown in Fig. 6. In the lower three pan- 
els, the dotted line shows the quasi-stationary X-ray component, 
and the dashcd-dottcd lines show the absorbed gamma-ray energy 
spectra normalized at 2 TeV to the observed ones. 



variable hard component are able to account for the large 

spectral changes observed at X-rays. 

Over the narrow spectral range form 3 keV to 25 keV 
we describe the quasi-steady X-ray component by a power 
law. We determined possible values of flux level and spectral 
slope of the quEisi-stationary component from extrapolating 
the 10 keV vs. 2 TeV flux correlation toward zero gamma- 
ray flux, and the 10 keV flux vs. 3-25 keV photon index 
correlation (Paper I) toward zero X-ray flux, respectively. 
Due to the scatter of both correlations this criterion gave a 
range of allowed values. We chose the values which resulted 
in the the best two-component SSC fits to the data, namely a 
10 keV amplitude v = 10~^° erg cm~^ s~^ and a photon 
index of 2.2. 

In the following we describe different incarnations of 
two-component models: two with variability through a time 
dependent rate of accelerated particles, and one with a time 
dependent Doppler factor of the SSC emission region. 

4-2. 1 Time variability through Qo{t), 1000 

First we consider a two-component model with flares caused 
by varying Q{){t). Since the spectral variability at X-rays 
is produced by the varying dominance of the soft quasi- 
static and the hard time dependent components no addi- 
tional spectral variability has to be produced by a changing 



7max and we use a fixed 7max corresponding to a high en- 
ergy cut-off in the synchrotron spectrum at MeV energies. 
Fig. 6 shows the fit of the two-component model (see fig- 
ure caption for model parameters). Compared to the one- 
component Qo{t) model shown in Fig. 2, the additional 
quasi-stationary soft component significantly improves the 
fit of the X-ray photon indices. The improvement is most 
pronounced for the April data. The model also describes 
well the TeV gamma-ray fluxes, with the notorious excep- 
tion of MJD 50544. 

Finally, we compare the model SEDs with the one mea- 
sured by BeppoSAX. Based on the data of one of the three 
BeppoSAX observations (we used the observation of April 
11), one can determine the spectrum of the quasi-stationary 
X-ray component outside the energy range covered by the 
RXTE observations. The other two observations can then 
be used to check the model predictions. The upper panel of 
Fig. 7 shows the electron energy spectra averaged over the 
integration time of the 3 BeppoSAX pointings. The lower 
three panels compare the modeled with the observed SEDs. 
By construction, the model describes the X-ray data of April 
7; the fit to April 11 is also good, but the model spectrum 
of April 16 is too soft. More detailed inspection shows that 
the model fails to describe the temporal evolution of the ^ 1 
keV fluxes, i.e. for April 16 it produces too much flux below 
1 keV. 
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Figure 8. Same as in Fig. 6, but for the two-component model 
with 7inin = 10^. The model parameters are: 5j = 45, R = 4.5 x 

10" cm, B = 1.1 G, tesc = 10-* iic-l, 7n,in = lO^ , Tmax = 

1.4 X 10'^, i = 0.5, n = 0.00. 

Remarkably, the predicted TeV energy spectrum, mod- 
ified by extragalactic extinction according to the LSDM 
model of Primack et al. 2001 fits the HEGRA time aver- 
aged spectrum very well. 



4.2.2 Time variability through Qo(t), "/„ 



A similar model with a high value of 7min ~ 10^ does not 
show these difficulties. In this case the break in the energy 
spectrum is more abrupt and the peak of the synchrotron 
SED of the time variable component is narrower than in the 
previous case. Viable models with high minimum Lorentz 
factors are located in a completely different region of pa- 
rameter space: at Doppler factor 45, we infer a magnetic 
field of B = 1.1 G and a radius of 4.5x 10^^ cm compared 
to the values of B = 0.014 G and R = 3 x 10^^ cm for 
the previous model. A large magnetic field is required to 
assure sufficiently rapid cooling of electrons with Lorentz 
factors above 7min to Lorentz factors below 7min. The latter 
electrons arc needed to produce the optical and UV seed 
photons, and partially, also for producing IC gamma-rays in 
the 250 GeV to ~ 1 TeV energy range. A small radius R fol- 
lows than from the requirement to produce the observed IC 
luminosity, given the large magnetic field and the "narrow" 
synchrotron SED. 

The best-fit result is shown in Figs. 8 and 9. In the lat- 
ter, it can be recognized that the use of = 10® substan- 
tially improves the fit to the broadband BeppoSAX data. 
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Figure 9. Same as in Fig. 7 but for the two-component model 
with 7inin = 10®; shown in Fig. 8. 



Even without any extragalactic extinction the model of the 
TeV gamma-ray data is very soft and only barely consistent 
with the observed data below 10 TeV. Only above 10 TeV, 
the model implies a slight amount of extinction. Note the 
pronounced break of the IC spectrum at ~2 TeV. Obviously, 
fitting a power law to a small portion of such a spectrum 
and inferring the DEBRA intensity from the deviation of 
the observed spectrum from this power law will not produce 
correct results. 



4-2.3 Variability through 5 X 10® 

In the framework of one-component models, an emitting 
blob with constant and isotropic emission in its rest frame 
but with a varying angle between its motion and the line of 
sight can not account for the 1997 X-ray and TeV gamma- 
ray fiares. The reason is that the large variability of the peak 
energy of the synchrotron SED would imply a large change 
of the blob's Doppler factor and as a consequence a much 
larger than observed fiux variability (Paper I). In a two- 
component model however, a variable Doppler factor can 
explain the flux variability: the X-ray energy spectra mainly 
change due to the relative dominance of the quasi-stationary 
and the time- variable X-ray components. Figs. 10 and 11 
show the two-component fit to the time resolved data and 
the broadband spectral data, respectively. The model gives 
an excellent fit to the data. This model is qualitatively very 
different from the other ones: time variability can be pro- 
duced on small time scales by changing Sj , and the electron 
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Date [ MJD-50000 ] 

Figure 10. Fit of a SSC model with two emission components: 
(i) a quasi-stationary X-ray component, and, (ii) a time variable 
X-ray/TeV gamma-ray component, flares produced through <5j(t). 
The upper panel shows here 5j (t) , and the other data and units are 

the same as in Fig. 2. The model parameters arc; R = 10^^ cm, 
B =0.16G,te3C =10i?C-l,7min = 4.5 X 10^ , 7max = 2.9xl0'^, 

^ = 0.5, T) = 0.04. 



spectrum is a steady state electron spectrum, and does not 
develop in time. We used here a rather low value of the 

Dopplor-factor, Sj = 10. As a consequence, after correction 
for DEBRA extinction, the TeV energy spectra are steeper 
than the observed ones. 



5 DISCUSSION 

In this paper we describe the time resolved modeling of 
the X-Ray and TeV gamma-ray data of a 2 month obser- 
vation campaign. The time resolved analysis is plagued by 
the sparse observational sampling and the unknown modi- 
fication of the TeV gammar-ray energy spectrum by extra- 
galactic extinction. However, modeling the X-ray fluxes and 
energy spectra and the relative changes of the TeV gamma- 
ray fluxes and photon indices allows us to exclude some hy- 
pothesis about the flare origin. Furthermore, we are able to 
verify that simple but self-consistently evolved SSC models 
based on canonical power-law energy spectra of accelerated 
electrons are able to account for the very detailed observa- 
tional data. More specifically, our conclusions from the time 
dependent modeling are as follows: 

(i) One-component models do not fully describe the data. 
While, by construction, the models succeed in accounting 
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/ |\ 
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Figure 11. Same as in Fig. 7 but for the two-component model 
with time variable <5j (t) shown in Fig. 10. 



for the temporal evolution of the X-ray fluxes they do not 
adequately predict at the same time the range of observed X- 
ray spectral indices, the broadband 0.1 keV-200 keV energy 
spectra, and the variation of the TeV gamma-ray fluxes. 

(ii) Two-component models give surprisingly good fits to 
the data. In these models, the X-rays originate from a super- 
position of a soft quasi-steady component and a hard rapidly 
variable component. We found two models which give an ex- 
cellent fit: in the first model flares are produced by a time 
dependent rate of accelerated particles. In the second model, 
a changing Doppler factor causes the flares. In both models, 
changes of the observed X-ray energy spectrum mainly re- 
sult from the relative dominance of the quasi-stationary and 
the time- variable X-ray component. 

(iii) Accurate fits to the BeppoSAX broadband data re- 
quire a large minimum Lorentz factor of accelerated particles 
on the order of 7min = 10^ . 

(iv) Degeneracy in both, model variant and jet param- 
eters, prevents us to use the time resolved SSC calcula- 
tions to substantially tighten the constraints on the amount 
of extragalactic TeV gamma-ray extinction, compared to 
earlier work (see e.g. Paper I, Guy et al. 2000, Vassiliev 
2000, de Jager & Stecker 2002, and references therein). The 
gammarray SEDs of Figs. 7 and 11 are consistent with the 
LCDM DEBRA model of Primack et al. (2001). In contrast, 
the model of Fig. 9, implies negligible extinction below ~ 
10 TeV. Especially the model with flux variability through 
the Doppler factor (Fig. 11) can produce very different in- 
trinsic gammar-ray SEDs while perfectly fitting the X-ray 
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1 10 

Enersv F TeV 1 

Figure 12. Change of GeV/TeV gamma-ray optical deptli as 
function of gamma-ray energy, inferred from comparing tlie model 
of Fig. 11 (lower solid lines, respectively) with the 1997 time av- 
eraged gamma-ray energy spectrum measured by HEGRA (Aha- 
ronian et al. 1999b, 2001a). The upper solid line is an upper 
limit for parameter combinations for which the 0.5-16 TeV en- 
ergy spectrum approximates a power law of photon index 1.5. 
The upper limits include 2cr statistical errors and take the 15% 
uncertainty of the absolute HEGRA energy calibration into ac- 
count. Allowed DEBRA models have to lie between the two solid 
lines. However, even if an absorption model lies between the two 
solid curves it does not imply that a valid SSC model exists such 
that the absorbed gamma-ray spectrum describes the data (the 
solid lines give a necessary but not sufficient condition that a 
DEBRA model is consistent with the SSC models and the data). 
The dashed and dotted lines show model predictions of Primack 
et al. 2001 ("LCDM" with Kennicutt and Salpeter Stellar Initial 
Mass Functions, respectively). All curves have been normalized to 
1 TeV where the systematic errors on the HEGRA energy spec- 
trum are small. 



and gamma-ray flux variations and the X-ray photon in- 
dices. In this model, the data constrain only the absolute 
flux level and energy spectrum of the quasi-stationary com- 
ponent, and the relative changes of the Doppler factor. The 
absolute value of , as well as the parameters R, B, and tesc 
remain degenerate. 

An upper limit on the modification of the TeV gamma-ray 
energy spectrum by extragalactic extinction can be derived 
from the fact, that the emitted time averaged gamma-ray 
energy spectrum is unlikely to be harder than dN^/dE cx 
with r fa 1.5. In Fig. 12 the range of allowed changes of 
gamma-ray optical depth with gamma-ray energy is shown 
and is compared to recent model calculations of Primack et 
al. (2001). A more accurate estimate of the amount of extra- 
galactic extinction from SSC modeling of Mrk 501 requires 
to pin down the jet parameters. Some key-observations are 
discussed further below; a more detailed discussion will be 
given by Coppi et al. (2002). 



(v) Table 2 lists for all studied models the electron to 
magnetic field energy density ratio r = {ue / Mb) as well as 
the minimum kinetic luminosity Lk = i^ cV^ (we -f wb) 
(Begelman et al. 1994) and we use F = All models are 
strongly out of oquipartition with r between 300 and 7500. 
Similar results, derived from an one-zone stationary SSC 
model, have recently been reported by Kino et al. (2002). 
The kinetic luminosities lie between 5 x lO''^ erg s~^ and 
2 X 10''* erg s^^, about 1000 times and more than the co- 
moving radiative luminosities which arc on the order of --^ 
5 X 10"^^ erg s~'. We computed the kinematic luminosity 
assuming a steady state jet with the same physical param- 
eters as the SSC emission region. Although this assumption 
might overestimate the true kinematic luminosity by a fac- 
tor of a few, our models clearly indicate that TeV blazars 
have rather powerful jets. Models with high 7min-values are 
closest to equipartition and require the leeist power. 

In SSC models the X-ray to TeV gamma-ray luminos- 
ity ratio strongly depends on the size and magnetic field of 
the emission region. As a consequence, most models that 
have been proposed to account for the flaring activity (as 
e.g. the internal shock model of Spada et al. 2001) do not 
naturally predict such a tight correlation of X-ray and TeV 
gamma^ray fluxes through a large number of distinct flares 
as evident in Fig. 1. The hypothesis that a single emission 
region of constant size produces a series of flares encoun- 
ters several problems: (i) due to the strong dominance of 
particle pressure over magnetic field pressure, the emission 
region should quickly expand adiabatically and thus become 
undetectable; (ii) it is not clear how the energy required for 
sustaining a prolonged flaring phase could be fed into the 
emission region; (iii) during the flaring period that lasted 
more than At ~ 2 months, the emission region would have 
advanced by ~ c F^ Af, that means by a distance of about 
~100 pc. The stability of the radius of the emission region 
would thus imply a jet opening angle of ~ 10"'* rad, several 
orders smaller than radio observations indicate. 

Our preferred interpretation is that flares originate from 
distinct emission regions with very similar characteristics, 
i.e. size and magnetic field. Such emission regions might form 
as the jet becomes radiative at a certain characteristic dis- 
tance from the central engine. The fact that our models give 
particle escape times on the order of and shorter than the 
flux variability time scale indicates that the flare duration is 
limited by the adiabatic expansion of individual emission re- 
gions. The jet would naturally feed energy to the site where 
the flares originate. The conclusions presented here are not 
limited to SSC models. Also in External Compton models, 
the tight X-ray /TeV gamma-ray correlation indicates a pre- 
ferred location for the production of individual flares: why 
else should the ratio of the jet frame magnetic field and ex- 
ternal seed photon energy densities remain roughly constant 
during 2 months? 

The preferred distance from the central engine could 
correspond to a characteristic length at which the jet be- 
comes unstable. Alternatively, a change in ambient pressure 
could induce jet instabilities at a characteristic distance from 
the central engine. Note that a qualitatively different but 
similarly puzzling stability has been found in the hardness 
intensity correlation of Mrk 421 (Fossati et al. 2000) for mea- 
surements taken between days and years apart. One conclu- 
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sion from this discussion is that refined modeling should 
treat adiabatic expansion in more detail. 

Since the modeling is computationally very intensive, 
we explored only a limited number of models. We did not 
consider External Compton models which historically have 
been applied to the more powerful EGRET blazars. For high 
jet Doppler factors even a very weak external photon field 
as e.g., IR radiation from dust, can be boosted and become 
significant in jet the frame. Depending on the seed photon 
energy spectrum, the radiative IC cooling of lower energy 
electrons might be stronger than for high energy electrons 
due to the Klein-Nishina effect. A possible consequence is 
that the energy spectra of External Compton models can be 
harder than for SSC models for a given value of B. Thus, 
the Klein-Nishina effect introduces a very rich behavior of 
External Compton models and the consequences of radiative 
cooling in the extremely "blue" TeV gamma- ray blazars can 
substantially differ from those in EGRET GeV blazars. 

Crucial advances in fixing model parameters will only 
be possible by substantially extending the observational cov- 
erage in time and wavelength. The 1997 April and May ob- 
servations had diurnal integration rates of typically 2 times 
20 min. Pinning down the evolution of the source during sev- 
eral fiares requires quasi-continuous monitoring over many 
days. Unfortunately, no sensitive X-ray all sky monitor with 
broadband spectroscopic capabilities will be available for 
the next several years or even longer. Such an instrument 
would be able to participate in intensive Multiwavelength 
campaigns on a large number of objects. The upcoming gen- 
eration of Cherenkov telescopes CANGAROO III, H.E.S.S., 
MAGIC, and VERITAS will have energy thresholds of be- 
tween 10 GeV and 50 GeV and a one order of magnitude 
higher sensitivity than present day instruments. The lower 
energy threshold is of crucial importance as it makes it pos- 
sible to asses the IC component at low energies where ex- 
tragalactic extinction is negligible (2 < 0.1) or much loss 
{z = 0.5 — 1) than at ~ 500 GeV. The new experiments 
should be able to reliably assess changes of the diurnal 
GeV/TeV energy spectra with a statistical and systematic 
accuracy in photon index of 0.05 or better, due to better 
gamma-ray statistics and improved detector calibration and 
atmospheric monitoring. Thus spectral changes as shown in 
Figs. 3, 4, and 10 will become measurable. 

A key observation for fixing the jet parameters is to 
measure a time lag between the X-ray and the TeV gamma- 
ray flux variability. A general prediction of SSC models is 
a time delay of approximately a light crossing time R 
between the leading X-ray and following gamma-ray fluxes. 
The measurement of this delay would allow one to deter- 
mine the size of the emission region. The requirement that 
the DEBRA reduces the 2 TeV flux by a factor of 5 or less 
would break the degeneracy in 5j and B. If the X-ray/TeV 
gamma-ray lag remains elusive it may bo that the determi- 
nation of the jet parameters of the TeV blazars detected so 
far has to wait until more reliable DEBRA estimates will be 
available derived from multiple blazar detections at redshifts 
between 0.05 and 1. 
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Table 1. Selected Blazar SSC Models Relevant to This Paper 



Authors 


Objects Studied 


Time 

Dependent? 


BED Peak Determined By 


Flare Mechanism 


DEBRA 
Extinction? 


Inoue & Takahara (1996) 


3C 279, Mrk 421 


No 


Cooling vs. Particle Escape 


Not specified 


No 


Bcdnarek & Protheroe (1997; 


Mrk 421, Mrk 501 


No 


Not specified 


Not specified 


Yes 


1999) 












Bottcher et al. (1997) 


Mrk 421 


No 




B, 7min 


No 


Mastichiadis & Kirk (1997) 


Mrk 421 


Yes 


Cooling vs. Particle Escape 


Qo, 7max, B 


No 


Plan et al. (1997) 


Mrk 501 


No 


7min 


7min5 7niax 


No 


Dermer et al. (1998) 


generic 


Yes 


Cooling vs. Particle Escape 


5j 


No 








and Plasmon Deceleration 






Chiaberge & Ghisellini (1999) 


generic 


Yes 


Cooling vs. Particle Escape 


Qo 


No 


Coppi & Aharonian (1999) 


generic 


Yes 


Cooling vs. Particle Escape 


Qo, B, 7max 


Yes 


Kirk & Mastichiadis (1999) 


generic 


Yes 


Cooling vs. Particle Escape 


Qo 


No 


Kataoka et al. (2000) 


PKS 2155-304 


Yes 


Cooling vs. Particle Escape 


'7m ax 


No 


Petry et al. (2000) 


Mrk 501 


No 


Cooling vs. Injection Time 


p 


No 








Scale 






Kusunose et al. (2000) 


generic 


Yes 


Cooling vs. Particle Escape 


■^m-a-u- ftilrrmiErli f.^^r> anH f-ar.^ i 
/max \_iji.ii.v^L*g,i.i ^esc "■'"'■'-i t'acc j 


No 


Tavecchio et al. (2001) 


Mrk 501 


No 


Not specified 


Change of 7b 


No 


Krawczynski et al. (2001) 


Mrk 421 


Yes 


Cooling vs. Particle Escape 


'Ymax 


No 


Sikora et al. (2001) 


3C 279, PKS 1406- 


Yes 


Cooling vs. Injection Time 


Qo 


No 




076 




Scale or 7inin 






Kino et al. (2002) 


Mrk 421, Mrk 501, 


No 


Cooling vs. Particle Escape 


Not specified 


No 




PKS 2155-304 










This work 


Mrk 501 


Yes 


Cooling vs. Particle Escape 


Qo, 7max, l5j 


Yes 








or 7i„in 







Table 2. Parameters of Models Shown in Figures 



Time Dependent 
Parameter 


Comments 


Si 


R [cm] 


BIG] 


tesc 


Tmin 


Tmax 




V 


Ue/UB 


[ergs-i] 


Qo{t) 


1-component 


45 


l.lxlQi^ 


0.014 


10 


1000 


2.5x10'' 


0.5 


0.2 


660 


1.2 xlO"-* 


7max(t) 


1-component 


45 


1.5x10^® 


0.0089 


3 


1000 


1.6-25x10® 


0.5 


0.2 


1200 


1.7x10-''' 


7max(t) <X Qo(t)^ 


1-componcnt 


45 


3.2xl0i5 


0.035 


3 


1000 


1.6-6.3x10*' 


0.5 


0.1 


860 


8.3xl0''3 


Qoit) 


2-component 


45 


3.4x101'^ 


0.014 


3 


1000 


2.3x10'' 


0.5 


0.4 


7470 


1.3x10'" 


Qo{t) 


2-component, 


45 


4.5xl0i3 


1.12 


10000 


1.0x10^ 


1.4x10'' 


0.5 


0.00 


290 


5.6 xlO-'^ 




high 7min 
























2-component, 


10 


lOiB 


0.16 


10 


4.5 xlO^ 


2.9 xlO''' 


0.5 


0.04 


2970 


2.8 xlO''^ 




high 7niin 
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